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The existence in the physical QCD vacunm of nonzero gluon condensates, such as {g^F'^), re¬ 
quires dominance of gluon fields with finite mean action density. This naturally allows any real 
number value for the unit “topological charge” q characterising the fields approximating the gluon 
confignrations which should dominate the QCD partition function. If q is an irrational number 
then the critical values of the 8 parameter for which CP is spontaneously broken are dense in R, 
which provides for a mechanism of resolving the strong CP problem simultaneously with a correct 
implementation of ?7a(1) symmetry. We present an explicit realisation of this mechanism within 
a QCD motivated domain model. Some model independent arguments are given that suggest the 
relevance of this mechanism also to genuine QCD. 
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I. INTRODUCTION 

The functional space Ta of gluon fields A with finite 
classical action (in the infinite volume limit), 

lim S'y[A] < oo, (1) 

y—»-oo 

can be divided into equivalence classes according to inte¬ 
ger values of topological charge 

^ = ( 2 ) 

Integer values originate from the purely topological prop¬ 
erties of the gauge group: the above integral can be 
rewritten as a surface integral and, since at infinity all 
fields can be at most pure gauge configurations, the value 
of that integral is determined by the gauge group alone. 
According to these equivalence classes, the QCD parti¬ 
tion function can be written in the form of a Fourier 
series 

CXD 

Z{9)= (3) 

V——00 

where is given by a functional integral over the helds 
belonging to the v-th sector of Ta- The representation 
Eq.Q is characteristic of non-abelian gauge groups and 
appeals to the existence of non-perturbative gluon field 
conhgurations, namely multi-instantons. As a matter of 
fact, from the very beginning the above construction is 
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based on a quasi-classical principle: the functional in¬ 
tegral is assumed to be dominated by fields with mini¬ 
mal classical action, while fields with infinite action are 
excluded from consideration ad hoc. However, as is ul¬ 
timately required by hadron phenomenology, the physi¬ 
cal QCD vacuum is characterised by finite nonzero gluon 
condensates Q, among which the lowest order conden¬ 
sate can be related to a mean action density times the 
coupling constant squared [2, 

{g^F;,,{x)F;,{x)) = 4 hm {g^Sv)IV ^ 0, (4) 

which can be provided only by the class of conhgura¬ 
tions already excluded from Eq. 0 , the held conhgura¬ 
tions with extensive (scaling with the volume V) classical 
action. One faces the necessity of dehning the QCD par¬ 
tition function as an integral over the space Fa of gluon 
conhgurations A which are allowed to have nonzero clas¬ 
sical action density at space-time inhnity. For instance a 
suitable requirement for Fa would be 

lim E“,,(a:)E;^(x) = const, (5) 

»-00 

with a nonzero constant on the RHS, which can be related 
to the trace anomaly of the QCD energy-momentum ten¬ 
sor ill Condition Eq. © produces no essential diffi¬ 
culties since the inhnity arising from the volume of the 
system is simply a matter of normalisation of the func¬ 
tional integral, or put differently, of normalisation of the 
vacuum energy. The actual values of the gluon conden¬ 
sates and the constant on the RHS of Eq. © are dehned 
by the minima of the quantum effective action. This in¬ 
dicates that the properties of the QCD physical vacuum, 
which are encoded in the various condensates and are 
responsible for conhnement and the mode of realisation 
of chiral symmetry, are due to purely quantum effects. 
Consequently, the helds to be identihed as dominating 
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the QCD functional integral must not have vanishing ac¬ 
tion density at infinity, otherwise the main effect, the 
existence of condensates, would be omitted. In order to 
define the functional integral over fields subject to Eq. 0 
in an analytical approach one can represent A = B + A 
with the background B from the class of fields dominat¬ 
ing the integral and A being small localised fluctuations 
in this background. The integral over fields B has to be 
performed exactly, while fluctuations A can be treated 
perturbatively. The candidate dominating fields B^{x) 
could be required to satisfy another condition for almost 
all X G R"' 

= const, (6) 

which, and this is important, neither forbids space-time 
variation of the strength tensor and the fields B^{x) 
nor devalues the importance of topological (singular) 
pure gauge field defects of various dimensions distributed 
throughout the entire Euclidean space-time Hi 

One can of course look for conditions other than 
Eqs. © and l©l to attempt to specify background fields 
B such that they would not be assumed to fill space-time 
almost everywhere and would be described in the spirit 
of dilute instanton gas or liquid models [^. However, de¬ 
spite diluteness, any superposition of an infinite number 
(in the infinite volume) of instantons and anti-instantons 
represents a field with infinite classical action, which can 
be different from the class fixed by Eq. ®, but for sure is 
absent in the space of integration subject to Eq. O- 
Eq. ® has an important consequence: in the space of 
fields with extensive classical action the functional 
Q{A\ can take any real value, rational or irrational, finite 
or infinite in the limit E —> oo. This signals that the 
equivalence classes fed into the representation Eq.® do 
not exist in the space Ta- Nevertheless, the notion of 
the mean value of the modulus |(5[-4]| 

Q = lim (|Qv[^]|), (7) 

V —^OO 

is well-defined on Ta and, as is known from phenomenol¬ 
ogy and lattice QCD, is a nonzero constant, related to 
the topological susceptibility y of the QCD vacuum 0. 
In both cases of fields satisfying Eq. ® and dilute su¬ 
perpositions of instantons, Q can be expressed in terms 
of a mean action density. There is a crucial difference 
however between instanton based approximations of the 
fields with infinite action and backgrounds (as in the do¬ 
main model for instance 0, IE S HEl ) based on Eq. ®. 
The literal use of instantons to approximate fields with 
infinite action assigns integer values of topological charge 
to the fields from Ta and thus transports integer units 
V of topological charge from Eq. ® into QCD with con¬ 
densates where such integers are irrelevant. Instead, ap¬ 
proximations based on Eq. ® leave the freedom to have 
any real values for Q for genuine gluon configurations and 
unit ’’topological charge” q for the fields approximating 
them in a particular model of QCD vacuum. 


It was realised a long time ago (see El , and, for a re¬ 
cent review ^El) that the strong CP-problem is a prob¬ 
lem of the quasi-classical treatment of the QCD func¬ 
tional integral built into the representation Eq. ®, and 
that incorporation of the long-range gluon configurations 
responsible for confinement could potentially remove this 
problem. However to date no way has been found to 
achieve this and maintain simultaneously the resolution 
of the Ca( 1) problem traditionally attributed to instan¬ 
tons. Scenarios based on the existence of an additional 
boson, the axiom have been suggested and have become 
canonical l[TE.[T^l. 

Nevertheless, as will be discussed below in detail on 
the basis of a QCD motivated model, and as has been 
proposed in 0 within the general model independent 
consideration, the strong CP and t/A(l) problems have 
a mechanism for simultaneous resolution within QCD. 
As is stressed in 0, the order of thermodynamic limit 
and the limit 0 —> 0 are in general not interchangeable 
and thus independence with respect to 9 of the infinite 
volume QCD partition function of QCD does not au¬ 
tomatically lead to vanishing topological succeptibility 
X = lime^o limy^oo V~^<PZv{9)/d9‘^. 

The QCD motivated model we shall consider in this 
article is the domain model 0 |E IE IS ■ For complete¬ 
ness and clear identification of the origin of the effect 
discussed in this article a comment about the choice of 
boundary conditions in the domain model is in order. 
The field B is considered as a background field and is 
subject to the condition Eq. ®. The system is consid¬ 
ered in a large but finite volume V. The total volume 
is split into a large number of subvolumes - domains. 
The fluctuation quark fields are defined by the bag-like 
conditions on the domain boundaries. The characteristic 
size of a domain is much smaller than the characteris¬ 
tic size of the total volume. The thermodynamic limit 
is understood as the limit when the size of the total vol¬ 
ume goes to infinity together with the number of domains 
while the size of domains stays finite. This construction 
is sufficient for defining and for practical calculation of 
the Euclidean functional integral for the partition func¬ 
tion. It should be stressed that bag-like conditions are 
very different from the (quasi-)periodic conditions typical 
for formulations based on a compactification of to the 
torus. The most important feature for our considerations 
is that the spectrum A of the Dirac operator under bag¬ 
like boundary conditions is asymmetric with respect to 
A ^ —A, while (quasi-)periodic boundary conditions lead 
to a symmetric spectrum. In particular, this asymmetry 
is responsible for the formation of the quark condensate 
and resolution of the Ua{F} problem in the model as is 
discussed in detail in 0. Nevertheless, there is no contra¬ 
diction between our results and canonical considerations 
based on torus-like compactifications of R^ M- These 
are just two different complementary statements of the 
problem. 

Within the domain model the mechanism for resolving 
the strong CP-problem is realised in the following way. 
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In the presence of nonzero gluon condensates irrational 
values of Q are permitted and lead to a realisation of CP 
in which the set of critical values of the 6 parameter for 
which the CP-breaking is spontaneous M is dense in the 
interval [—7r,7r]. As a consequence, the infinite volume 
partition function with infinitesimally small quark masses 


Z = lim Zv{d) = lim Zv{0) 

V —>-oo V—»-oo 

is independent of 0, and 


hm (E)t. = hm (E)^=o, 

V —»-oo V —»-oo 

lim = 0, 

V —>-oo 

for any CP-even and CP-odd operators E and O respec¬ 
tively, which resolves the problem of CP-violation. Si¬ 
multaneously one finds that topological succeptibility in 
QCD with massive quarks 


is nonzero, independent of 0, and satisfies the anomalous 
Ward identity, which indicates a correct implementation 
of the Ua{^) symmetry. In particular, in the chiral limit 
the mass of the rj' is expressed via the topological suc¬ 
ceptibility of pure gluodynamics in agreement with 
the Witten-Veneziano formula. 

If Q takes a rational value then only a finite number 
of such critical points exist in the interval 9 € [— 7r,7r], 
which results in the standard strong CP problem. 

To apply a mechanical analogy, an irrational value of Q 
leads to a picture which is reminiscent of ergodic motion 
resulting in a dense winding of a torus by a trajectory 
|T^. while motion over closed trajectories can be associ¬ 
ated with rational Q. 

This article is devoted to an investigation of this mech¬ 
anism on the basis of the domain model of QCD vac¬ 
uum iiilli We omit here all details related to 
the motivation, formulation and numerous results of the 
model given in the above papers, apart from those which 
are specifically needed here. In the next two sections the 
9 dependence in the model under consideration is eluci¬ 
dated. The Witten-Veneziano mass formula, Gell-Mann- 
Oakes-Renner relation, the anomalous Ward identity and 
r/TTTT decay are specifically discussed in section ITvl The 
final section is devoted to explanation of the relation be¬ 
tween the mechanism of the strong CP resolution in the 
domain model and the Poincare recurrence theorem of 
ergodic theory which justifies the title of the present pa¬ 
per. 


II. THETA DEPENDENCE IN THE DOMAIN 
MODEL 


The model is given in terms of the following partition 
function for N oo domains of radius R 

Z = Af lim f dadp 

V,N^od / 

X J )A«]AfpA«] 

where the functional spaces of integration and are 
specified by the boundary conditions {x — ziY = B? 

(x) = 0, 

i fji{x)e''^°‘'^^ ^ ='ip^''\x), ( 9 ) 

Here fii = with the generators of 5't/c(3) in the 
adjoint representation, the a and [3°“ are random chi¬ 
ral angles associated with the chiral symmetry violating 
boundary condition Eq. m in the presence of Nf {a = 
1,..., Nj — 1) quark flavours. The A“ are SU (Nf) flavour 
generators normalised such that Tr(A“A^) = 25°"^. Aver¬ 
aging over the angles a, /3“ is included in the partition 
function. The thermodynamic limit assumes V, N oo 
but with the density v~^ = N/V taken fixed and finite. 
The partition function is formulated in a background field 
gauge with respect to the domain mean field, which is ap¬ 
proximated inside and on the boundaries of the domains 
by a covariantly constant (anti-)self-dual gluon field with 
the field-strength tensor of the form 

N 

t=i 

with = B'^Si.p. Here z^ are the positions of the 

centres of domains in Euclidean space. The measure of 
integration dai over parameters characterisin g d omains 
can be found in any of the papers 11 ISIS The 
background field in Eq. is designed to satisfy Eq. ©■ 
In a sense, this model represents a “step-function” ap¬ 
proximation to a class of fields with infinite action which 
are assumed to dominate the QCD partition function. 
The building blocks of this approximation are domains 
with the same mean size R and mean action density B^, 
whose values are determined from the phenomenological 
string tension in . Moreover as the field in domains is 
taken to be (anti)-self-dual, a mean absolute value q of 
the integral of the topological charge density over the do¬ 
main volume V = 7r^i?'^/2 is attributed to a domain, and 
is expressed through the domain radius and field strength 
a.s q = B^R'^/16. This quantity can take any real val¬ 
ues but for B, R fixed as in g « 0.15 .... It should 


f[J dai J 
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be noted also that domain boundaries are assumed to 
be populated by pure gauge singularities which look like 
topological defects (instantons, monopoles, vortices and 
domain walls) but only locally; globally no topologically 
conserved number can be associated with them. The role 
of these defects in the model is to generate boundary con¬ 
ditions for the fluctuations of gluon and quark fields A 
and Ip inside domains 0 . In particular, at the location of 
pure gauge singularities dVi colour quark currents have 
to satisfy the condition 

=0, xG dVi. 

with rj'p being a vector normal to the surface dVi of the 
j—th domain where the singularities are located. This 
leads to the condition Eq. 0 with arbitrary parame¬ 
ters a, /3“. Simultaneously it is natural to require that 
colourless quark currents like 

'ip{x)X°^’tp{x),'ip{x)j 5 ^‘'P^{x),'ip{x) fji{x)^ 5 X°‘ip{x), etc, 

are continuous at dVi. Such continuity conditions 
Ea. dlOII fix the parameters a, /3“ to be the same at all 
singular surfaces. 

Apart from chiral angles and neglecting gluon fluctu¬ 
ations A in lowest order of the perturbation expansion, 
the integrations in the partition function Eq. © give 0 

Zv{0) = N I dad0exp{—VJ^(a,/3,d)}, 

flclS 

with the free energy density as a function of chiral angles 
in the presence of infinitesmally small quark masses rrii 

j Nf 

T = -ln[cOS(7(ITAr, — 0)1 — TOiCOS$i, (10) 

i—1 

and where 

Nf 

arctan(tan<i)i), 4)^ = a + Bi. 

i=l 

The functions Bi for various numbers of flavours are 

Si = 0 for Nf = 1, (11) 

Si = B 2 = for Nf = 2, (12) 

and, rather than give explicit representations for Nf = 3, 
we give several identities useful for later purposes: 

3 

y^s, = 0 ( 13 ) 

\ (14) 

2=1 a —1 

= i4E(/3“)'/^®+--- (15) 

2=1 ^ a —1 


where the dots denote other cubic terms but involving 
a = 4,...,8. Ens. dTar^ come from the diagonalisation 
of the fermionic boundary conditions Eq. 0 in the flavour 
space. The quantity H has been explicitly computed in 
and appears in the domain model as a function of the field 
strength S and domain radius R 

1 °° h 

where M is the confluent hypergeometric function, z = 
fiBR^ / 2 and the colour Tr denotes summation over ele¬ 
ments of diagonal color matrix n. With S, R as deter¬ 
mined in i H = (237.8MeV)^. It determines the value 
of the quark condensate 

{'ipp;) = —H 

and arises from the asymmetry of the Dirac operator 
spectrum with the boundary conditions Eq. @ . Due 
to Ea. lfl^ (which applies for any Nf) we have 

Wnj = Nfa (mod tt). 

Thus the anomalous part of free energy density, that in¬ 
volving the topological charge g, is independent of the 
/3“ which manifests the Abelian property of the anomaly. 
The significance of Eq. is that it is derived directly 
from the domain ansatz for the QCD vacuum but also 
corresponds to the zero momentum limit of the effective 
chiral Lagrangian of 0 up to the somewhat different en¬ 
coding of the anomaly. A task in the following is then 
verifying that this encoding of the symmetries of low en¬ 
ergy phenomenology reproduces the expected phenomena 
of spontaneous chiral symmetry breaking and resolution 
of the Ua(X) problem but more importantly that it also 
solves the strong CP problem within this QCD vacuum 
approach. 

In the thermodynamic limit E ^ oo the minima of 
T define the vacua of the system under consideration. 
For massless quarks rrii = 0 the anomalous term gives a 
discrete set of degenerate minima at 

= + + ^ {k&Z,leZ). ( 16 ) 

Thus due to the anomaly the continuous Ua{^) is reduced 
to a discrete group of chiral transformations connecting 
the minima Eq. C0 with each other. In the absence of 
additional sources violating chiral symmetries like quark 
masses, this discrete symmetry is sufficient to provide for 
zero mean values of all chirally non-invariant operators. 

III. IRRATIONAL VALUE OF PARAMETER q 

There is a crucial difference between rational and irra¬ 
tional values of q. For q rational there is a finite number 
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of distinct minima and all others are 27r-equivalent to 
one of these minima. In other words, there is a periodic 
structure which allows one to split the set Eq. into 
finite number of equivalence classes. 

If q is an irrational number then all minima in the in¬ 
terval a € [—(X),oo] are distinct: there are no minima 
which are 27r-equivalent to each other. The distance be¬ 
tween a pair of minima is 


A(fci, ^ 2 ) = 2Trm{ki,li\k2,l2) + S{ki,h\k2,l2) 


with m{ki,li\k 2 ,l 2 ) some integer, and ^(fci, ^i|fc 2 , ^ 2 ) a 
number in the interval ( 0 , 27r) which is different for any 
pair (fci,/i|fc 2 , ^ 2 ), h ^ h- Unlike the case of ratio¬ 
nal q where only several numbers S arise, the set of all 
i5(fci, Zi|fc 2 , Z 2 ) is dense in the interval ( 0 , 27r) for irrational 
q. The periodicity, typical for a rational < 7 , is lost for ir¬ 
rational values of q and the set of all minima cannot be 
split into 27r-equivalence classes. However, these discrete 
minima of the free energy are separated by infinite en¬ 
ergy barriers in the infinite volume limit, and no flavour 
singlet Goldstone modes are expected. 

For any real q this discrete symmetry and contin¬ 
uous non-singlet flavour chiral symmetry are sponta¬ 
neously broken, which becomes manifest if the infinites¬ 
imal masses are switched on. The term linear in masses 
in Eq. selects a minimum for which aki{0) and the 
flavour non-singlet angles /3“ maximise ^ cos(<I>i). 
There are specific values of the 6 which are critical in the 
sense that two different minima aki{dc) and ak'i'{0c) are 
degenerate in the presence of a mass term; these degen¬ 
erate minima are CP-conjugates. The so-called Dashen 
phenomenon occurs [13: CP is broken spontaneously for 
0 = 9c while for non-critical values CP is explicitly bro¬ 
ken. For rational q there is only a finite number of 9c 
in the interval (—tt, tt], which includes HHl the value 
9 = TT. Irrational values of q lead to a drastically dif¬ 
ferent picture: the set of critical values 9c is dense in K. 
In other words, any real value of 9 is either critical or 
is a limit of a Cauchy sequence of critical values. CP 
is broken spontaneously in any arbitrarily small vicin¬ 
ity of any real value of 0 , in particular in the vicinity of 
9 = 0. Understood in the sense of limiting values, the 
pseudoscalar condensate vanishes in both CP-conjugate 
vacua. The qualitative (and, in a sense, dramatic) differ¬ 
ence in the consequences of q being rational or irrational 
is illustrated in FigUl 

Note, that at 0 = 0 CP is just an exact symmetry, it is 
not broken at all (neither spontaneously, nor explicitly). 
The point 0 = 0 is a limiting point of a set of critical 
points but it does not belong to this set: 0 = 0 itself 
is not critical. Thus there is no contradiction with the 
Vafa-Witten theorem [T^. 





0 0 


FIG. 1: The scalar (A and B) and pseudoscalar (C and D) 
quark condensates as functions of 9 for Nf = 3 in units of H. 
The plots A and C are for rational q = 0.15, while B and D 
correspond to any irrational q, for instance to q = ^ = 

0.15047.. . which is numerically only slightly different from 
0.15. The dashed lines in A and C correspond to discrete 
minima of the free energy density which are degenerate for 
m = 0. The solid bold lines denote the minimum which is 
chosen by an infinitesimally small mass term for a given 6. 
Points on the solid line in A, where two dashed lines cross each 
other, correspond to critical values of 9, at which CP is broken 
spontaneously. This is signalled by the discontinuity in the 
pseudoscalar condensate in C. For irrational q, as illustrated 
in B and D, the dashed lines densely cover the strip between 1 
and —1, and the set of critical values of 9 turns out to be dense 
in R. As will be discussed in the last section, an analogy can 
be seen between A (C) and motion on a torus over a closed 
trajectory, while a manifestation of the Poincare recurrence 
theorem resulting in dense winding of a torus (for example, 
see fl3 l can be recognised in B (D). 

IV. MASS RELATIONS AND ANOMALOUS 
WARD IDENTITY 

Flavour SUh{Nf) x SU^iiNf) is broken spontaneously 
and the appropriate number of corresponding Coldstone 
modes are expected, while Ua(1 ) is realised in a more 
complex non-Coldstone form for any 0, as shown in M- 
To make this manifest at the level of the free energy 
Eq. (1101 one can consider small variations of the variables 
a and in the vicinity of the minimising configurations 
{a = aniini9),Bi = B 2 = = 0}. To this end it is 

sufficient to consider variations 5a, 5(3°". We organise the 
expansion as follows: 

T = Tfj + Ta- ( 17 ) 

The leading term is Tq = — N f^m coslaki (9)] which 
we have already identified in as generating the 
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quark condensate per flavour including its 9 dependence, 
{QnQn) = —^cos[afei(0)]. Fixing now Nf = 3 we next 
consider the terms involving fluctuations in (5/3“: 


Tp = -Hmcos[afci(6>)] ^(5/3' 


a\2 


a—1 


--^^'^msm[aki{6)]^{6(3‘^fS/3^ (18) 

^ a=l 


where we have showed all quadratic terms and one of the 
cubic terms. Eqs. ma are used to derive this result. 
Let us now identify the flavour structure of the fermionic 
boundary condition Eq. m with that used in effective chi¬ 
ral theories 0 via 

z(/3“A72 + a) ^ 

where 7 r“ and rjo are respectively identified as the non¬ 
singlet and singlet mesons and the t“ are normalised such 
that Tr(r“r“) = 6°^^. The pion-decay constant is = 
93 MeV. Thus we identify, for example for the neutral 
flavourless mesons, the fluctuations as 




27ro 

F ’ 

TT 


TT 


Sa = 



(19) 


Certainly what follows is not more than an illustration 
of what one can expect for the meson spectrum and in¬ 
teractions: a more solid investigation requires the con¬ 
struction of the effective meson action and calculation of 
the relativistic bound state spectrum in the spirit of 0 . 
Nevertheless, if we identify the coefficients of the terms in 
Eq. HI 811 which are quadratic in the variations as “meson” 
mass terms then we obtain 


= 2 mH. 

This takes into account that for irrational q one has 
cos[afei( 6 >)] = 1 and sin[Q!fci( 6 >)] = 0 for any 9 (see the 
bold solid lines in plots B and D in FigP). Using our 
identification of H with the quark condensate, we thus re¬ 
cover the Gell-Mann-Oakes-Renner relation for masses 
of pions (a = 1, 2, 3) and the eta-meson (a = 8 ). 

More significantly for the strong CP-problem we ob¬ 
serve that in higher order terms in the decomposition 
those which are CP-even, namely involving even powers 
of the fluctuations, come with a coefficient cos[afe;(d)] = 1 
for the allowed vacua. All CP-odd terms, those with odd 
powers of the fluctuations, are proportional to sin[afe;(d)] 
which exactly vanishes. CP is thus unbroken. For exam¬ 
ple, the term which would correspond to an rjinr inter¬ 
action is the cubic term in Eq. m- Using the identifica¬ 
tions Eq. m, this is 


- sin[afc;( 0 )] 


2V3F, 


-TTsiTTof. 


( 20 ) 


For q irrational this term exactly vanishes: the decay 
r] ^ inr is exactly suppressed. For q rational, we can 
expand Eq. (HU in the vicinity of d = 0 using Eq. 
and obtain — /Ft^)tts{ttq)'^ which is the analogue 

of Eq. ( 8 ) in (and similar results elsewhere) for three 
de gene rate quarks. In the context of approaches such 
as |2flj, there is an a priori assumption that one can ex¬ 
pand in small 9: the 9 dependence is first chirally rotated 
from the topological charge density term into the quark 
mass terms and then an expansion performed for small 
quark mass. Using the small upper bound of the neu¬ 
tron dipole moment such considerations lead to the 
expectation of an unusually small value for 9, leading to 
the familiar fine-tuning problem. However, in the context 
of irrational topological charge in the domain model, as 
a particular study of the class of fields specified by the 
condition Eq. m where the 9 dependence can be handled 
exactly, we see that such an expansion misses the possibil¬ 
ity of the cancellation of CP-conjugate vacua in CP-odd 
quantities. The absence of CP-violating hadronic pro¬ 
cesses such as 77 —> TTTT (or, if baryons could be included 
in these considerations, a neutron dipole moment) is not 
then indicative of vanishing 0 ; physical observables are 
independent of 9, as in Figs^B,D), and still CP symme¬ 
try is exact in QCD with condensates, confinement and 
the correct realisation of chiral symmetry, even in the 
presence of the 9 term. 

Einally we turn to the fluctuations in the singlet direc¬ 
tions a. Consider then the integral 

r+00 

Zv{9)=N / da 


with V = vN and finite ^ 1. The integral can be 
computed by means of the saddle-point approximation. 
The set of minima of the free energy density is obtained 
from the solution of equation 

J-'^(a, Pinin,9) = —v~^ ta,ii(qNfa — 2qkTr — q9 — 2lTr) 

—Nfm't^ sm{a) = 0 . 

The free energy density and its second derivative at the 
minima aki is 

F = —'AmNfCos[aki{9)] 

Nj+ '^mNfCOs[aki{9)] 

where = (f'/'v = / 123-n'^ is the topological 

susceptibility in the absence of quarks which in the do¬ 
main model was evaluated in 0 to be approximately 
(197 MeV)"*. With the identification of 6a of Eg. 11911 
we thus extract the Witten-Veneziano mass formula for 
the eta-prime meson 

ml,Fl = 2Nfx'^^+mlF^^ (21) 

in the physical vacuum. Continuing with the evaluation 
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of the integral in the standard way we arrive at 


V. DISCUSSION 


Zv{0) = lim 

L—>oo 

^^ -A^^l ^NNfvm'H, cos[a.ki{0)] 

S + HmNfCOs[akii0)]]^/^ 

X [l + C>(TO)][l + C>(l/iV)], (22) 


with ZviO) = 1- The sum above runs over all local min¬ 
ima of the free energy, enumerated by k, 1. 

Recall now that the topological susceptibility 


X = 


9 


327r2 J 


d^x{F ■ F{x)F ■ F(0)) 


= — lim lim 


1 1 


Zv{0)V 


Zv{0). 


(23) 


It can thus be extracted from the partition function by 
the double variation with respect to 9. It is important to 
note here that this relies on changing the order of taking 
the derivatives with respect to 9 and functional integra¬ 
tion, which is well-defined only for a certain regularized 
form of the functional integral. In particular the total 
volume V must be large but finite. In our context it 
means that the differentiation of Eq. (1221 with respect to 
9 must be performed before taking the thermodynamic 
limit iV —> oo. When V is large but finite Zv{9) de¬ 
pends on theta. After the second derivative over 9 is 
taken the limit N —>■ oo chooses the global minimum, 
i.e. l,k corresponding to the global minimum. The key 
point here is that for irrational g, V e > 0 and V 9 there 
exist integers {l,k) (specifying the vacuum aki(9)) such 
that I — cos{9/Nf + 2'Kl/q/Nf) < e. The existence of 
such a minimum eliminates dependence of observables 
on 9. Moreover, taking the limit d —> 0 is unnecessary in 
Eq. (EHll : as defined by Eq. 123) but with d ^ 0 omitted, 
the CP-even quantity y does not depend on 9. We re¬ 
peat though that the most important point here is that 
differentiation with respect to 9 and the thermodynamic 
limit are not interchangable; but deriving with respect to 
9 first is precisely how y is obtained through Eq. (EHll . 

Now let us use this to extract the topological suscepti¬ 
bility in the presence of light quarks. The aim, on the one 
hand, is to check consistency with the anomalous Ward 
identity of , and on the other to confirm the suppres¬ 
sion of the susceptibility by quark loops, as proposed by 
123 to fulfill the identities. According to Eq. ll^ the 
result is : 


= ^ + (24) 

and consistent (in Euclidean space) with the anomalous 
Ward identity of [^: 

Nj^QCD ^ NfmlFl + 0{ml). (25) 


A. The Poincare recurrence theorem 

It is instructive to compare some aspects of the pic¬ 
tures arising for rational (in particular, integer) versus ir¬ 
rational values of q. The fundamental difference between 
these two cases is clearly seen in Fig. ^ but still requires 
precise formulation. As a matter of fact, plots (A,B) rep¬ 
resent all inequivalent curves — cos[Q;fc;(^)] while (C,D) 
show all inequivalent curves sin[afci( 0 )]. 

First of all, for rational q there is a finite number 
ric of critical points ({ 6 >c}, i = I,...,nc) in the inter¬ 
val —TT < 9 < n and a countable set of critical points 
({0* -I- 27rj}, i = l,...,nc,j G Z) in R, distributed 
(27r/nc)-periodically in R. In particular, 0 = 0 is not 
critical while d = tt is critical for any rational q. The 
number of critical values ric is determined by the integer 
numerator qi of a rational number Nfq = qxjq^- For 
irrational q the interval —tt < 9 < already contains a 
countable set of critical points ({ 0 ^}) * = 1, • ■ •, oo) which 
is dense in this interval. Thus V0 G R we have either a 
critical point or a limiting point of a sequence of critical 
points. In particular, d = 0 is such a limiting point for 
any irrational q. 

Furthermore, for rational q and V0, the set {Cjifi)} 
(j = 1,..., ric) contains a finite number of points. Here 
Cj{9) = cos[aj(0)]. The representation for Zv{9) is 
analogous to Eq. (EH) but contains a finite sum over 
j = I,..., ric) only one term of which dominates the sum 
V 0 ^ 9c\ the rest of the terms are exponentially sup¬ 
pressed at large volume. Such a representation indicates 
that Zy(9) is differentiable to infinite order V 6 > ^ 9c^ 
which in particular means that CP-even observables can 
be computed for small 9 as Taylor series in 9. At critical 
9 the first derivative of Zv{9) is discontinuous. 

Irrational values of q change this behaviour radically: 
C = I is the limiting point of a sequence {Cj{9)} {j = 
1,..., oo) V0. This means by definition that Ve > 0 and 
V0 3Je(6() such that 1 — Cj{9) < e, Vj > Je{9). Simulta¬ 
neously |5'j(6()| < e, Vj > Je{9). Here Sj{9) = sin[aj( 6 ()]. 
The dependence of Jg on 9 is discontinuous V0 and thus 
the sum in Eq. 1221 ) in the infinite volume limit represents 
the function Zv{9) which, strictly speaking, is not differ¬ 
entiable V0. However, for finite volume, each term in the 
series is a differentiable function and the sum, regularised 
by L < oo and N < oo, can be differentiated term by 
term, with a subsequent taking of the limits V, L —>■ oo 
which are well defined. In this sense Eq. ( 1221 ) gives a 
prescription for computing derivatives with respect to 9. 

To underscore that the character of the 9 dependence 
realised in the domain model does not represent anything 
unusual or “exotic” it is useful to map the above consid¬ 
erations to some textbook examples of application of the 
Poincare recurrence theorem - one of the corner stones 
of ergodic theory (for example, see p^l. 

The first straightforward example is just the present 
problem slightly rephrased. Consider the group of rota- 






FIG. 2: Values of 9i = g^9 as defined in Eq. 12611 : boxes occur 
for rational q = qi/q 2 with qiNf = 3, while circles represent 
some of the points densely distribnted, in particular, in the 
neighborhood of 6^ = 0 for irrational q. 

tions by angle f3 acting on a point at polar angle 0 on a 
circle. The Z-th element of this group corresponds to the 
angle 

9i = g^9 = 9 + (31 (mod 27r), (3 = 2^9. 

According to the Poincare recurrence theorem, for any 
rational b there exists a number Zf, such that 0;^ = 9. 
The points represented in Fig.[51by the boxes correspond 
to the elements of the group. For any irrational b the 
set {9i, I = 0.. .00} is dense on the circle everywhere, 
particularly in the vicinity of 0 = 0 as is illustrated by 
the filled circles in Fig. |21 Identifying the 6*-parameter of 
QCD with the polar angle in Fig. [3 and the parameter b 
with the “topological charge” qhy b = 1/Nfq establishes 
the desired map. 



FIG. 3: Coordinates on a torus, as discussed in the text. 

Another simple example of the application of the 
Poincare theorem is that of uniform motion on a torus, 

0 = </* = b^, 

where 9 and (j) are the latitude and longitude of a point P 
on the torus as is shown in Fig.|3 If the ratio of velocities 
b = be/b^ is a rational number then the point moves 
along a closed trajectory which can be characterised by 
an integer winding number. In the case of irrational b the 
trajectory is dense on the torus. If we take for simplicity 
6 ^ = 1, identify bg = 3/qNf and replace the discrete 
parameter 2ttI by a continuous variable t, then the plot 
A in Fig^would correspond to motion on the torus over 
a closed trajectory with a winding number equal to the 
number of critical points ric- An irrational q leads to a 


correspondence between the curves densely covering the 
strip in plot B to the motion of a point on a torus over 
the trajectory densely covering the torus. 

If we were to take the liberty of extrapolating the QCD 
motivated domain model considerations of this paper to 
genuine QCD, then we would conclude that for the class 
of background fields with Q\A\ taking rational values the 
Ua{3) problem can be solved but at the cost of simulta¬ 
neously inviting the strong CP problem. On the other 
hand, irrational values of Q[A] are as natural for QCD 
with nonzero gluon condensates as rational (in particular, 
integer) values but providing clear signatures of resolving 
the Ua{3) problem without creating the CP problem in 
strong interactions. 

B. Model independent consideration 

In this final part of the discussion we give arguments 
as to how the mechanism for resolving the strong CP- 
problem can be identified in QCD in a model indepen¬ 
dent way. We start with the anomalous Ward identities 
of, for example, M- These identities are model inde¬ 
pendent and arise from the symmetries of QCD and ba¬ 
sic properties of the pseudoscalar meson spectrum. A 
compact summary of this approach is given in Appendix 
B.l of M- Without repeating those derivations here, 
a key result obtained is a differential equation for the 9 
dependence of the matrix of condensates 

Vij oc {0\qiqj\0) (26) 

whose phases are parametrised by angles (pi. The true 
vacua are obtained by finding local minima with respect 
to perturbation of these angles by small 

From this one derives that a CP-even symmetry break¬ 
ing term is given by 

{(^Kven) = ^ COS (p, (27) 

i 

and a CP-odd symmetry breaking term is given by 

= 2 to ^ sin . (28) 

i 

One has the following model independent requirements 
then: 

No Z7(l) Goldstone boson, summarised in the equation 

i 

Stationarity condition with respect to chiral perturba¬ 
tions, reflected by the equation 

sin(()i = independent of i. (30) 

Minimum with respect to chiral perturbations, 
cos (pi > 0. 


(31) 
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CP- conservation, 

^ sin (j)i = 0, (32) 

i 

which might hold for only one value of 6 (eg 0 = 0 ) or, 
ideally, may somehow hold for all 9. 

2^-periodicity. This is summarised in the equation 

^ cos (()i I e = ^ cos I e+ 2 n 7 i- (33) 

i i 

but cosine could just as well be replaced here by sine. 
This condition emerges from previous conditions but only 
in the presence 0 of the Dashen phenomenon at least 
at 0 = TT. 

Each of these conditions is realised in the domain 
model ^3 has appeared in previous sections of this 
paper. In the following we shall only derive results from 
these equations. 

From Eq. we extract the result 
(/), = 0 + c 

i 

for some 0-independent constant c. We can give the in¬ 
dividual angles as 

(j)i = {9 + Ci)/N, N~^ Cj = c. 

i 

From Eq. (EOl we have in turn that 

sin[(0-|-ci)/A^] = sin[(0-|-C2)/A^] = • • • = sin[(0-|-CAr)/iV] 

thus 

Ci = C = KTT, 

is independent of i up to shifts of 2k7T. All the allowed 
solutions to Eqs. i29l3l)l31ll are classified by 

4>i = ''' = 4 *n = + ‘^Trk 

with the restriction ipo = {ktt 0)/N e] — 7 r/ 2 , 7 r/ 2 [, 
providing for positivity of the cosine. At this point k is 
a real number. 

Let us first implement CP-invariance at 0 = 0 just 
keeping k arbitrary. Since (pi are independent of i it 
suffices that sin^j = 0. Thus 

sin KTT/N = 0 

which leads to k = Nm. with m integer. Thus either k 
can be absorbed in k or should be set to zero. The latter 


choice leads to the statement in 0 that CP-invariance 
at 0 = 0 fixes unambiguously the integration constant 
from Eq. ra . 

Alternately let us characterise k before seeking to im¬ 
pose CP-invariance at 0 = 0. Note to this end that any 
real number can be written as a product of an integer 
and some number between zero and one, 

k = 2/^, ZgZ, Ce [0,1/2]. (34) 

This opens the possibility that different values of I are 
appropriate for different ranges of 0 subject to the above 
conditions. Thus far for CP-invariance at 0 = 0 we must 
have 1 = 0. But for 0 7 ^ 0 nonzero values of I satisfying 
condition 

21^71 -1-0 = 0 (mod 27r). 

are now available due to the discrete nature of this label. 

We are now ready to map the model independent ap¬ 
proach identically to the framework that emerged from 
the domain model via the the presence of the parameter 
which can have irrational values. In view of Eq. (Unj it 
is tempting to make the identification 

? = 1/q. 

The conclusion we draw is that the spectrum of solu¬ 
tions to the theta dependence of the QCD Ward identi¬ 
ties is rich enough to exhibit the behaviour seen in and 
originally derived in an explicit QCD motivated model, 
the domain model, which solves the strong CP-problem. 
However the interpretation of the parameter is possible 
at this stage only within the model, namely as the mean 
value of a fraction of the topological charge related to the 
space-time regions where the dominating gluon configu¬ 
rations in the QCD vacuum can be considered as homo¬ 
geneous, or succinctly - the mean topological charge per 
domain. 
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